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We consider the properties of a gas of bosonic diatomic molecules in the limit when few of the 
molecules are dissociated. Taking into account the effects of dissociation and scattering among 
molecules and atoms, we calculate the dispersion relation for a molecule, and the thermal depletion 
of the condensate. We calculate the dependence of the Bose-Einstein condensation temperature of a 
uniform gas on the atom-atom scattering length, and conclude that, for a broad Feshbach resonance, 
the condensation temperature increases as the molecular state becomes less strongly bound, thereby 
giving rise to a maximum in the transition temperature in the BEC-BCS crossover. We also argue 
on general grounds that, for a gas in a harmonic trap and for a narrow Feshbach resonance, the 
condensation temperature will decrease with increasing scattering length. 
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I. INTRODUCTION 



The crossover between BCS pairing of two species of 
fermion and Bose-Einstein condensation of a diatomic 
molecules with changing atom-atom scattering length has 
been studied by many authors, beginning with the works 
of Eagles |jy] , Leggett ^ and Nozieres and Schmitt-Rink 
One of the surprising results of the calculations of 
Ref. H and also found in Refs. [J, |5j is that, in the 
crossover region, the transition temperature has a maxi- 
mum as the strength of the atom-atom interaction is var- 
ied. However, it is unclear whether the maximum is real, 
or merely an unphysical artifact of either the assump- 
tions made in the model or the numerical procedures. 
The early calculations only accounted for quadratic fluc- 
tuations around the mean-field solution and it was gen- 
erally expected that the fluctuation induced peak was an 
artifact of the gaussian approximation. This view was 
furthermore supported by a self-consistent field theoret- 
ical calculation in Ref. 'u] where Tc was found to be a 
monotonic function of the dimensionless coupling con- 
stant. However, the most recent self-consistent calcula- 
tions of Tc do indicate the presence of a peak in Tc- 
Most calculations for broad resonances for a uniform gas 
give a peak in Tc for both one-channel and two-channel 
models for the interaction, and for a variety of approx- 
imations for the atom propagators, which are in some 
cases replaced by either bare propagators, or by dressed 
propagators 0, la> 1^ ■ In calculations that use a combi- 
nation of bare and dressed propagators, Tc initially de- 
creases coming from the far BEG limit before turning 
around to give a local maximum near resonance (IL] . In 
calculations for a harmonically trapped gas in the local 
density approximation, the peak vanishes for both one- 



channel and two-channel models for the interaction, and 
for a variety of approximations for the atom propagators 
[1^ . The peak vanishes in calculations for sufficiently 
narrow resonances. A summary of the present situation is 
that microscopic calculations point to there being a max- 
imum in Tc for broad resonances for the uniform gas, but 
not for harmonically trapped gases or sufficiently narrow 
resonances. 



The purpose of this paper is to study the properties of 
a weakly dissociated gas of diatomic molecules composed 
of fermionic atoms with a view to determining whether or 
not the transition temperature Tc to a paired state should 
exhibit a maximum in the BEC-BCS crossover. We shall 
consider the molecule effective mass and the reduction of 
the density of molecules due to dissociation, as well as 
the interaction between molecules. As a consequence of 
critical fluctuations, molecule- molecule interactions lead 
to an increase in Tc away from the BEC limit and this ef- 
fect dominates the decreases in Tc we find due to changes 
in the effective mass and the number of molecules. For a 
trapped gas, we shall argue that Tc should decrease away 
from the BEC limit due to the reduction of the central 
density in the trap as a consequence of the increasingly 
repulsive molecule-molecule interaction. 



This paper is organized as follows. In Sec. II we con- 
sider the uniform system, and calculate the effective mass 
of a molecule and thermal depletion of the number of 
molecules. The transition temperature is discussed in 
Sec. HI, first for a uniform gas and then for a trapped 
gas. Finally, Sec. IV contains concluding remarks. 
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II. THE UNIFORM SYSTEM 

We investigate the properties of a system containing 
equal numbers N of two species of fermion, labelled by 
a pseudospin index ct = ±1, that can form a diatomic 
molecular bound state in the limit when the system con- 
sists mainly of molecules, with a small admixture of 
atoms. We shall assume that the two species of fermion 
have the same mass m. In the spirit of Landau, we write 
the effective low-energy Hamiltonian for the system as 
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the operators a\ a , b create and destroy atoms and 
molecules, respectively. Here V is the volume of the sys- 
tem, g is the matrix element for decay of a molecule into 
two atoms and aij is the scattering length between a par- 
ticle i and a particle j, and 'a' stands for an atom and 
'm' for a molecule. The quantity aaa is the background 
scattering length for atom-atom scattering due to nonres- 
onant processes. In the case of strong coupling between 
the molecule and the atoms, e. g. for a broad Feshbach 
resonance, and in the limit of a deep ly bound molecular 
state it has been shown that 
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(2) 



We begin by considering a uniform gas, and we first 
calculate the molecule spectrum. The leading contribu- 
tions to the energy of a molecule are due to the molecule- 
molecule interactions, which give a shift 2Trh?nn^aynm/'n^ 
to the self energy of a molecule in the mean field of the 
other molecules, and a contribution one half of this to the 
average energy per molecule. Here is the number den- 
sity of molecules. Atom-molecule interactions contribute 
dirfi^niiaam/m to the molecule self energy from the scat- 
tering of atoms by molecules, and from dissociation of 
molecules into atoms an amount 



T-p/2+q 



^p/2-q 



{2TThf eb + (q + p/2)2/^" 



The term independent of the atom distribution function 
Uq is not included, since this is already contained in the 



binding energy Cb of two fermions in the absence of a 
medium. The contribution ||2J) is due to Pauli blocking 
of dissociation of molecules, which raises the energy of 
a molecule. As a consequence of Galilean invariance, in 
the absence of a medium, the effective mass of a molecule 
remains equal to the bare mass. The binding energy of 
a molecule with momentum zero is thus reduced by an 
amount 



{2TThY eb + q^/m 



25^-, 

eb 



(4) 



where ria is the density of atoms of both species. In the 
second expression we have assumed that the molecule is 
deeply bound in the sense that eb ^ kT. The molecule 
spectrum varies over momentum scales of order (meb)^^^, 
and therefore for the momenta of interest in the prob- 
lem of Bose-Einstein condensation, which are of order 
(mfcr)^/^, it is sufhcient to calculate only to order p^. 
The result is 
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(5) 

where the effective mass of a molecule is increased by an 
amount given by 
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(6) 



We now turn to the transition temperature of the gas. 



III. TRANSITION TEMPERATURE 

A. Uniform system 

For a uniform gas of elementary bosons, the transition 
temperature is unaffected by the Hartree mean-field en- 
ergy, since this gives a shift to the energy of a boson and 
to the chemical potential which are independent of the 
momentum of the particle. However, it has been demon- 
strated that, due to critical fluctuations, the transition 
temperature is increased, by an amount which is given 
by IIllll 



for small rim'^amrn 



5Tc 1 /o 
— = bnj^ an 



(7) 



where 6 is a positive coefficient, which has been esti- 
mated to be approximately 1.3 0, 0, |23|. This ef- 
fect will tend to be counteracted by the increase of 
the effective mass and the reduction of the density of 
molecules due to dissociation, both of which tend to de- 
crease the transition temperature, since for an ideal Bose 
(3) gas, Tc ~ ft^nm'^/2M*. The number density of molecules 



is given by 



N 
V 



2 



(8) 
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The two latter effects give changes in the transition tem- 
perature which are proportional to n^juy^ oc 
Thus for a deeply bound molecule, the effects of crit- 
ical fluctuations dominate, and the transition tempera- 
ture initially increases away from the BEC limit. Since in 
the BCS limit the transition temperature tends to zero, 
we conclude that the transition temperature must have 
a local maximum, provided that the transition tempera- 
ture is a continuous function of the scattering length. 

B. Trapped cloud 

We now turn to the case of a cloud in a harmonic trap, 
and we shall assume that the trapping potentials for the 
two species of fermion are the same, and equal to one half 
of the trapping potential for a molecule. Previously, the 
transition temperature for a trapped gas of bosons has 
been calculated with allowance for the nonzero spacing 
of the single-particle levels in the oscillator potential and 
for the distortion of the density profile by the interaction 
and the result is |23| 

^ ^ _0.73— 7V-i/^ - 1.33^7Vi/6^ (9) 
Tc a 

where is the algebraic mean of the trap frequencies for 
the three principal axes of the trap, (D is their harmonic 
mean and a = (H/mujy^'^ is the mean harmonic oscillator 
length. 

To take account of critical fluctuations, we argue 
that provided the typical interaction energy between 
molecules, 2TrH'^ni^anun/'m' is large compared with the 
oscillator quantum of energy but small compared with 
kTc, it should be a good approximation to calculate the 
change in Tc in a local density approximation, in which 
the central density of the cloud is inserted in Eq. Q. 
The central density nni(O) at Tc is given by 

n^(0) = C(3/2) [p^J , (10) 

and since /cTc = [N/C{'i)Y/^hbj, it follows that the con- 
tribution to the change in Tc due to critical fluctuations 
is 

<^Tc 1.3 [C(3/2)]V3 i/ea^n, 

Tc ~ (2^)1/2 [C(3)]i/6 a ■ ^ ' 

Calculations indicate that critical fluctuations are sup- 
pressed in inhomogeneous systems |22| , so we regard this 
estimate as an upper limit to the increase of Tc. Adding 
this to the earlier result © without fluctuations, we find 
for the total change in Tc 

^ « _0.73— - 0.64^iVi/6. (12) 
Tc ui a 



Thus, in contrast to the uniform case, the change in Tc 
for small is negative, since the effect of critical fluc- 
tuations is overwhelmed by the reduction of the central 
density. 

The discussion we have given above applies to broad 
Feshbach resonances, for which the energy dependence 
of the atomic interactions may be neglected. For narrow 
Feshbach resonances, critical fluctuations are suppressed 
because the magnitude of the contribution from processes 
in which the argument of the interaction contains a fre- 
quency which is to be summed over will be reduced be- 
cause of the strong reduction of the interaction away from 
resonance, as discussed in Ref. |23l[ . 



IV. CONCLUSION 



In this paper we have identified a number of effects that 
produce changes in the transition temperature away from 
the BEC limit. Two of these are due to thermal disso- 
ciation of molecules, and these lead to a decrease in Tc. 
However, for a uniform gas, these effects are overwhelmed 
by a rise in Tc due to the increasingly repulsive molecule- 
molecule interaction. We therefore conclude, quite gen- 
erally, that Tc must display a global maximum between 
the BEC and BCS limits for a broad Feshbach resonance. 
For trapped gases, we have argued that Tc will initially 
decrease away from the BEC limit. Thus, on the basis of 
the arguments given here, there need not be a maximum 
in Tc but it cannot be excluded. Likewise, for a narrow 
Feshbach resonance, critical fluctuations are suppressed, 
and consequently Tc is expected to decrease away from 
the BEC limit. 

The framework we have described in this paper could 
be a useful tool for analysis of detailed crossover calcu- 
lations, since it points to a number of quantities that 
could be calculated in order to test our fundamental un- 
derstanding. 
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